nucl-th/0306033 



p meson photoproduction at low energies 

Yongseok Oh 1 ' 2 '* and T.-S. H. Lee 2 '* 

Institute of Physics and Applied Physics, 
Yonsei University, Seoul 120-749, Korea 
2 Physics Division, Argonne National Laboratory, Argonne, Illinois 60439 

(Dated: February 4, 2008) 

Abstract 

The o"-exchange and /2-exchange mechanisms for p meson photoproduction are re-examined. 
Then the commonly employed <7-exchange amplitude is revised by using the recent information 
from the analyses on the p — > 7r°7r°7 decay and the aNN coupling constant from Bonn potential. 
Instead of relying on the Pomeron-/ proportionality assumption, the meson exchange amplitude 
is established from an effective Lagrangian which is constructed from the tensor structure of the /2 
meson. Phenomenological information together with tensor meson dominance and vector meson 
dominance assumptions are used to estimate the f2 coupling constants. As a first step to improve 
' the current theoretical models, we have also explored the effects due to the un-correlated 2tt 

exchange amplitude with ttN intermediate state. This leading-order 2tt exchange amplitude can 
be calculated using the coupling constants determined from the study of pion photoproduction 
and the empirical width of p — > tttt decay. In comparing with the existing differential cross section 
data, we find that a model with the constructed 2ir, a, and ji exchanges is comparable to the 
\Q • commonly used a exchange model in which the a coupling parameters are simply adjusted to fit 

the experimental data. We suggest that experimental verifications of the predicted single and 
double spin asymmetries in the small \t\ (< 2 GeV 2 ) region will be useful for distinguishing the 
two models and improving our understanding of the non-resonant amplitude of p photoproduction. 
Possible further improvements of the model are discussed. 
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I. INTRODUCTION 



The recent experiments at Thomas Jefferson National Accelerator Facility (T JNAF) [1-4] , 
GRAAL of Grenoble [5], and LEPS of SPring-8 [6] are expected to provide new opportunities 
for studying the electromagnetic production of vector mesons at low energies. For example, 
the differential cross section data for p photoproduction from the CLAS Collaboration at 
TJNAF show big differences with the old data of 1970's [7, 8] in the large momentum- 
transfer (|t|) region at low energies, where one may learn about the VNN couplings and 
other production mechanisms [9-12]. Much more new data with similar high precisions will 
soon be available. 

The study of vector meson photoproduction is expected to shed light on the resolution of 
the so-called 'missing resonance' problem [13-17]. On the other hand, it is well known that 
this can be achieved only when the nonresonant mechanisms are well understood [18, 19]. As 
a continuation of our effort in this direction [14, 18], we explore in this work the nonresonant 
mechanisms of p photoproduction. 

There exist some investigations of the nonresonant mechanisms for vector meson pho- 
toproduction. To account for the diffractive features of the data in small t region at high 
energies, the Pomeron exchange model, as illustrated in Fig. 1(a), was developed. However, 
this model fails to describe the experimental observables at low energies. Indeed, meson 
exchanges (or secondary Reggeon exchanges) are found to be crucial in understanding the 
low energy data. In the case of uj photoproduction, it is well known that one-pion exchange 
is the most dominant process at low energies. For p photoproduction, however, the situa- 
tion is not clear. Generally, there are two scenarios which are based on either the a meson 
exchange model [20, 21] or the ji exchange model [10, 22]. The a exchange model was 
motivated [20] by the observation that the decay width of p — > ixix^f is much larger than 
the other radiative decays of the p meson. It is further assumed that the nn in the nnj 
channel can be modeled as a a meson so that the p<yy vertex can be defined and modeled for 
calculating the a exchange mechanism as illustrated in Fig. 1(b). In practice, the product of 
the coupling constants g pal gcrNN of this tree-diagram is adjusted to fit the cross section data 
of p photoproduction at low energies. If we use g%. NN /4ir ~ 8 from Bonn potential [23], we 
then find that the resulting g pai will yield a decay width of p° — > o"7 an order of magnitude 
larger than the value extracted from the recent experimental decay width of p° — > 7r°7r°7 
[24-26] . Thus the dynamical interpretation of the commonly used a meson exchange model 
for p meson photoproduction must be further examined theoretically. 

In this work, we would like to take a different approach to account for the exchange of 
ixix in p photoproduction. First, the commonly employed cr-exchange amplitude is revised 
by using the coupling constant g a NN from Bonn potential and g pal from the recent exper- 
imental decay width of p° — > 7r°7r°7 with the assumption [20] that 7r°7r° in this decay is 
strongly correlated and can be approximated as a a particle. This is our starting point of 
developing a new model which is more consistent with the existing meson-exchange models 
for NN scattering [23], irN scattering, and pion photoproduction [27]. We then consider the 
consequence of the strong p — > 7r + 7r~ decay which accounts for almost the entire decay width 
of the p meson. With the empirical value of the p meson decay width, one can define the 
pixix vertex, which then leads naturally to the 'uncorrected' two-pion exchange mechanism 
illustrated in Fig. 2 with M = n in the intermediate state. A more complete calculation 
of un-correlated 27r-exchange contributions to p photoproduction should also include other 
intermediate states such as uN and tcA. However, the contributions from these interme- 
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diate states involve propagation of two or three pions and must be considered along with 
other multi-pion exchange mechanisms (such as the crossed diagrams due to the interchange 
of 7 and p lines in Fig. 2). Obviously, this is a much more complex task and will not be 
attempted in this exploratory investigation. Our calculation of 27r-exchange will be detailed 
in Section II. F. 

The f'2 exchange model for p photoproduction was motivated by the results from the 
analyses of pp scattering data at low energies [28] . In the study of pp scattering the dominant 
secondary Regge trajectory is represented by the / trajectory, and the idea of Pomeron-/ 
proportionality had been used to model the Pomeron couplings using the ji couplings until 
1970's [29-32] before the advent of the soft Pomeron model by Donnachie and Landshoff 
[33]. By considering the role of the / trajectory in pp scattering, it is natural to consider 
the /2 exchange model for vector meson photoproduction. However, the ji exchange model 
developed in Refs. [10, 22] for p photoproduction made use of the Pomeron-/ proportionality 
in the reverse direction. Namely, they assume that the structure of the f?, couplings are the 
same as that of the soft Pomeron exchange model. Thus the fi tensor meson was treated as 
a C = +1 isoscalar photon, i.e., a vector particle. In addition, the fit to the data is achieved 
by introducing an additional adjustable parameter to control the strength of the f 2 coupling 
[10]. This is obviously not very satisfactory and leaves a room for improvement. 

Instead of relying on the Pomeron-/ proportionality assumption, the f'2 meson exchange 
amplitude is evaluated in this work starting with an effective Lagrangian which is con- 
structed from the tensor structure of the /2 meson. Phenomenological information together 
with tensor meson dominance and vector meson dominance assumptions are used to es- 
timate the f'2 coupling constants. With this, we then explore the extend to which the p 
photoproduction data can be described by a model that includes this newly constructed /2 
exchange amplitude together with the revised a-exchange amplitude and the un-correlated 
27r-exchange amplitudes discussed above. 

This paper is organized as follows. In Section II, we explicitly define the amplitudes for the 
considered p photoproduction mechanisms, including the Pomeron exchange, a exchange, 
pseudoscalar meson exchanges, s- and w-channel nucleon terms, and the newly constructed /2 
exchange. The 2-7T exchange amplitudes are then given to complete our model construction. 
The numerical results are presented in Section III. For comparison, we consider two models. 
Both models contain the s- and -u-channel nucleon terms and the exchanges of Pomeron, 
7r, and rj. In addition, the first model includes the a exchange with free parameters to fit 
the data following Refs. [20, 21], while the second model contains the two-pion, a, and 
/2 exchanges, where the parameters of the a exchange are fixed by Bonn potential and 
p — > 7r°7T 7 decay. We explore the extent to which these two rather different models can 
be distinguished by examining the differential cross sections and spin asymmetries. Section 
IV contains a summary and discussions. The details on the / 2 interactions with the photon 
and hadrons are given in Appendix for completeness. 

II. MODELS FOR p PHOTOPRODUCTION 

In this Section, we discuss possible production mechanisms for 7p — > pp. We first dis- 
cuss single particle exchanges as depicted in Fig. 1. Then the 2n exchange model will be 
constructed. Each of the considered production amplitude, as illustrated in Fig. 1, can be 
written as 

T fi = el{V)M^eM, (1) 
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FIG. 1: Models for p photoproduction. (a,b) i-channel Pomeron and one-meson exchanges (M = 
f2,7r,7],a), (c,d) s- and u-channel nucleon pole terms. 



where sJV) and £^(7) are the polarization vectors of the vector meson and the photon, 
respectively. We denote the four-momenta of the initial nucleon, final nucleon, incoming 
photon, and outgoing vector meson by p, p', k, and q, respectively. The Mandelstam variables 
are s = W 2 — (k + p) 2 , t — (p — p') 2 , and u — {p — q) 2 . 



A. Pomeron exchange 



We first consider the Pomeron exchange depicted in Fig. 1(a). In this process, the incom- 
ing photon first converts into a qq pair, which interacts with the nucleon by the Pomeron 
exchange before forming the outgoing vector meson. The quark-Pomeron vertex is obtained 
by the Pomeron-photon analogy [33], which treats the Pomeron as a C = +1 isoscalar pho- 
ton, as suggested by a study of nonperturbative two-gluon exchanges [34]. We then have 
[33, 35-37] 

= G P (s,t)T^, (2) 

with 

r r = ii Wy^ 1 ., Ui?%- t ) ww)it,r - w), (3) 

where a em = e 2 /An and Fi is the isoscalar electromagnetic form factor of the nucleon, 

F J t) = 4 < - 2M (4) 

1{) (4M2 r -t)(l-t/0.71) 2 ' {) 

with t in GeV 2 . The proton and vector meson masses are represented by M N and M v , 
respectively. (M v = M p in our case.) 
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The Regge propagator for the Pomeron in Eq. (2) reads 

(\ a p (t)— 1 r ■ \ 

-J expj-^MO-1])- (5) 

The Pomeron trajectory is taken to be the usual form a P (t) = 1.08 + a' P t with a' P = l/s = 
0.25 GeV~ 2 [33]. In Eq. (3), fy is the vector meson decay constant: f p = 5.33, f u = 15.2, and 
/ = 13.4. The coupling constants f3 u = (3 d = 2.07 GeV" 1 , (3 S = 1.60 GeV" 1 , and /i 2 = 1.1 
GeV 2 are chosen to reproduce the total cross section data at high energies, E 1 > 10 GeV, 
where the total cross section of vector meson photoproductions are completely dominated 
by the Pomeron exchange. For p photoproduction, we set (5 q = (3 q r — (3 U — (3d- 



B. a meson exchange 

The a meson exchange model advocated by Friman and Soyeur [20] is based on the 
observation that Y{p — > nit^) is the largest among all p meson radiative decays, which 
leads to the assumption that the p photoproduction process at low energies is dominated 
by the exchange of 2ir. The In is then effectively represented by a o meson. The effective 
Lagrangian for this model reads [9, 20, 21] 

C a = ^ (PfTdrA, - d»p»d u A,) a + g aNN NaN, (6) 
p 

where p M is the p° meson field and the photon field. The resulting a meson exchange 
amplitude is 



>C = 69p ^ NN - ^ m (k ■ qgT ~ u(p')u(p)F„ NN (t)FUt), (7) 

where 



A 2 a - Ml Ap CT7 - Ml 

FaNNif) = ^ 2 _ ^ , F pal {t) = — — (8) 

are the form factors. The cutoff parameters of the form factors and the product of coupling 
constants g pal goNN are adjusted to fit the p photoproduction data at low energies. It was 
found [20, 21] that 

M a = 0.5 GeV, gl NN /*n = 8.0, g pai = 3.0, 

K = 1.0 GeV, A pCT7 = 0.9 GeV. (9) 

The resulting a mass parameter is close to the value M a = 0.55 ~ 0.66 GeV of Bonn potential 
[23]. If we further take the value ^ATAf/^ 71 " — 8-3 ~ 10 from Bonn potential, we then find that 
the resulting g pal is close to the values from the QCD sum rules, g pcr7 (QCDSR) = 3.2±0.6 [38] 
or 2.2 ± 0.4 [39]. However such a large value of g pal corresponds to the p — > o"7(^ 7r°7r°7) 
decay width that is much larger than the empirical value of T(p° — > 7r°7r°7) [24, 26]. If 
we accept the empirically estimated, but model-dependent value of SND experiment [26], 

BR(p -> (77) = (1.9 -oi ±0.4) x 10~ 5 , which gives T{p -> (77) w 2.83 keV, we get 

|^| w 0.25, (10) 
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since the Lagrangian (6) gives 

T(P--7) = ^^(M P 2 -M CT 2 ) 3 . (11) 

This value is smaller than that of Eq. (9) by an order of magnitude. Therefore, the a 
exchange model suffers from the big uncertainty of g pa ^i which is under debate [24-26, 40, 41]. 
Furthermore, there is no clear particle identification of a a particle and the use of a exchange 
in defining NN potential has been seriously questioned. Thus it is possible that the a 
exchange may not be the right major mechanism for p photoproduction. 



C. Pseudoscalar meson exchanges 

The 7r and rj meson exchanges are also allowed for p photoproduction, although their 
contributions are known to be not important. They are calculated from 

My 

C*NN = fg^75^iV, (12) 

where ip = n°,r]. The coupling constants g pw are fixed by the p — > yry decay widths 

n P ^^) = ^f(M v -Mlf. (13) 

Using the experimental data [42], T(p° — > 7r°7) cxpt . = 121 ± 31 keV and r(p° — > ^7) cxp t. = 
62 ± 17 keV, we get 

g PJ7T = 0.756, 9(rm = 1.476. (14) 

This also gives <7 W77r = 1.843 and g^ ir] = 0.414. We use gl NN /4ir = 14.3 and the SU(3) 
relation to get g^ NN /4n = 0.99. Although there are other estimates on the value of g V NN 
reported in the literature, the role of the r\ exchange is much suppressed in p photoproduction 
and the dependence of our results on g n NN is negligible. 

The pseudoscalar meson exchange amplitude, Fig. 1(b), calculated from the Lagrangian 
(12) reads 



i e 9(ry<p9'pNN 1 

where the form factors are 



M * = 2M„My t-M* ^ qak ^U ~ t'Hu(p)KNN(t)F pvi (t), (15) 



A 2 -Ml A 2 - M- 

F<pNN{t) = ^ 2 _ f P , F fxpj(t) = ^2 _ * ■ (16) 

We use A,, = 0.6 GeV, A pnJ = 0.77 GeV, A,, = 1.0 GeV, and A W7 = 0.9 GeV [14, 20]. 
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D. Nucleon pole terms 



The s- and it-channel nucleon terms, Figs. l(c,d), are calculated from 



C 



ipp 



-ppp 



-eN 

9pNN 



N 



-^—a lxv d v A* i 
2M N 



N, 



2M 



N 



N. 



(17) 



The resulting production amplitude is 



M7 = 



M 



N 



+ r»(k) * J+% N r$(q)F N (u)]u(p), 



u-M 2 N 



where 



F M - ^ - ? p cr^n 



/Cr. 



2Mjv 

The form factor has the form [43] 



^ N + (r-M* N y 



(18) 
(19) 

(20) 



with A at = 0.5 GeV taken from Refs. [9, 14]. This choice of the nucleon form factor leads to 
a satisfactory explanation of the steep rise of the differential cross sections with increasing 
\t\ in terms of the it-channel nucleon term [Fig. 1(d)]. 

Because Fn(s) ^ Fn(u), the above amplitude does not satisfy the gauge invariance. In 
order to restore the gauge invariance, we project out the gauge non-invariant terms as 



r fi -pp. 



(21) 



For the pNN coupling constants, we take the values determined in the analyses of pion 
photoproduction and 7riV scattering [27]: 



9 p nn = 6.2, k p = 1.0, 
and the anomalous magnetic moment of the nucleon is k p — 1.79. 



(22) 



E. / 2 meson exchange 

We now discuss the exchange of the / 2 (1270) tensor meson, which has quantum numbers 
I G (J PC ) = 0+(2++). The mass and decay width of the / 2 (1270) are M f = 1275 A ± 1.2 

MeV and r(/ 2 ) = 185.1 ^6 MeV [42]. Because of its quantum numbers, it has been once 
suggested as a candidate for the Pomeron. But this assumption violates the duality with 
the a 2 trajectory which includes I G (J PC ) = 1~(2 ++ ) state and it is now believed that the 
/ 2 does not lie on the Pomeron trajectory. 
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G 2 fNN /^ 


FfNN/G fNN 




1.12 


— 


Ref. [44] 


3.31 


W 


Ref. [47] 


3.31 ±0.63 


0.06 ±0.17 


Ref. [48] 


4.0 ±1.0 


0.00 ± 0.07 


Ref. [49] 


2.2 ±0.9 


0.6 ±0.9 


Ref. [50] 


0.38 ±0.04 


« 


Ref. [46] 



TABLE I: Estimates on the fNN coupling constants, GfNN and F/atat, using ttN dispersion 
relations. The values are compared with the prediction of tensor meson dominance [46]. 



In the approach of Ref. [10], the f 2 is treated as a C = ±1 isoscalar photon just like the 
Pomeron. This leads to a Regge amplitude of the following form 



Mf 2 =K f G f2 (8,tyip, (23) 



where 1 



'' W 2sin[xa /2 (t)]r[ a/l (t)] 

with s 1 — l/a'f 2 ~ 1 GeV 2 , while the form of T£ v is the same as given in Eq. (3). The f 2 
trajectory is linearly approximated as ctf 2 (t) « 0.47 + 0.89t [22, 28]. In order to control the 
strength of the f 2 couplings to the hadrons, a free parameter n f was introduced [10] and 
adjusted to fit the p photoproduction data at low energies. 

In this paper, we depart from this Regge parameterization and construct an f 2 exchange 
model solely based on the tensor structure of the f 2 meson. We will use the experimental data 
associated with the f 2 meson, the tensor meson dominance, and vector meson dominance 
assumptions to fix the f 2 coupling constants, such that the strength of the resulting f 2 
exchange amplitude is completely fixed in this investigation. Following Refs. [44, 45], the 
effective Lagrangian accounting for the tensor structure of the f 2 NN interaction is written 



as 2 



C fNN = -2i%™ JV( 7/ A + l„d»)Nr + 4^d,Nd„Nr, (25) 
where f^ u is the f 2 meson field. This gives the following form of the fNN vertex function, 

V be , d = -efufa) |^ [E m7 „ + 7 „E„] + ^f^} «(P6), (26) 

where S M = (p& + pa)^, Pb and pd are the incoming and outgoing nucleon momentum, 
respectively, and e^ v is the polarization tensor of the f 2 meson. 



1 The form of Gf 2 in Eq. (24) is due to the fact that the / 2 interaction is treated as that of an isoscalar 
photon, i.e., a vector particle interaction. If we use the tensor structure of the / 2 interaction, it would be 



«/ 2 (*)- 2 (l + cxp[-iwa f {t)])ira' 



h 



2sin[7ra /2 (i)]r[a /2 (t)-l]' 
In the conventions of Ref. [45], G^ NN — Gf NN and G^ N = Ff NN . 
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The coupling constants associated with the /2 meson were first estimated by using the 
dispersion relations to analyze the backward ttN scattering [44] and the tttt — > NN partial- 
wave amplitudes. The results are summarized in Table I. Note that the value estimated 
based on the tensor-meson dominance [46] is much smaller than the empirical values. (See 
Appendix for details.) 

The most general form for the /V7 vertex satisfying gauge invariance reads [51] 

( 7 (*MW 2 > = ^V^(M'), (27) 

where e and e' are the polarization vectors of the photon and the vector meson, respectively, 
and 

A%L(k, k>) = ^ [g KX (k ■ k') - k'M (k - k'),(k - k% 

+ 9fv y [9K\(k - fc%(fc - k') v + g Xp k' K (k - k') v + g\ v k' K {k - fc% 

- 9^k x {k - k') u - g KV k x (k - fc% - 2k ■ k'(g KfM g X u + fiWfiv)]- (28) 

The tensor meson dominance assumption together with the vector meson dominance gives 
[51] 

e 

f fVl = and g fVl = —G fvv , (29) 

where 

Gfvv = G f7T7T = 5.76. (30) 

Here Gf wn is determined from the decay width of ji — > tttt. The details on the interactions 
with the photon and hadrons, and tensor meson dominance are given in Appendix. 
With the above formulas, it is straightforward to obtain the production amplitude as 

M% = -u(p)r*W)ti(p) {p M f pa ^(k, q)F fNN (t)F fv ,(t), (31) 

where 

1 1_ 
Pa/3;pa ~ -{ga P g/3a + 9aa9/3p) ~ -^9a/39pa, 

V^(k,q) = lj^[-g^(k-q)+q^](k + q) p (k + q)„ 

- gvpk^k + q) a - g va k^(k + q) p + 2k ■ q(g up gpa + g va g^ , (32) 



+ 



9nu(k + q) p (k + q) a - g w> q v {k + q) a - g P *qu(k + q) f 



and 



9^-9 ^-^-^' . (33) 
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J(k) 



piq) 




M(q') 





-r 




! n 


Pip) 


B(p B ) pip') 



FIG. 2: 2ir exchange in p photoproduction. The intermediate meson state (M) includes tt, and the 
baryon (B) includes the nucleon. 



The form factors are chosen as 



fNN 



(*) = 



A 



fNN 



M 



h 



A 2 - t ' 

fNN 1 



Ff V7 it) - 



(34) 



where the cutoff parameters will be discussed in Sect. III. The relative phases among / 2 
couplings are fixed by tensor meson dominance. 



F. 27r exchange 

In this subsection, we discuss the 2-7T exchange for p photoproduction as shown in Fig. 2. 
We only consider the contribution from setting the intermediate state (MB) = (ttN). As 
discussed in Section I, the contributions from other intermediate states like nA and uN 
involve propagation of two or three pions and hence are neglected along with the other 
multi-meson exchange amplitudes in this exploratory investigation. 

We compute the loop amplitude of Fig. 2 by making use of the method of Sato and Lee 
[27], which gives 



Ti oop = J rf 3 q' [e( 7 ) • B jN , MN ik, q'; E)] GW(q, E) [e*(V) ■ V MN>pN (4, q : E)\ , 



(35) 



where 

Gmn(*{,E) = e _ E N {q') - E M (<fi + in ' ^ 

Obviously, 5 7 jv,mjv and Vmn, p n are the one-pion-exchange amplitudes illustrated in Fig. 2. 
We only consider (MiV) = (irN) intermediate state in this paper. 
Equation (35) can be rewritten as 



T loop = V I rf 3 q' 



J s*( V)-V(q!, tl )s^)-B(k,q') 
W - E B iq>) - E M (q') 

- i I dtt kt p BM (k t ) e*iV) ■ V(kt, q) e( 7 ) • B{k, k t ) 6{W - M M - Mb), (37) 
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where the subscripts of V and B are understood. Here 6(x) is the step function and 

where Es{k) and EM(k) are the energies of the intermediate baryon and meson with mo- 
mentum k. Through the on-shell condition W = E B (k t ) + E M (k t ), k t is determined as 

kt = ^KW 2 ,Ml,MH (39) 

where 

X(x, y, z) = x 2 + y 2 + z 2 - 2(xy + yz + zx). (40) 

For the considered (MN) = (irN) case, the one-pion-exchange amplitudes B^^n and 
V n N,pN m Eq.(35) can be calculated from 

X 7r|, A,, 

^iv = f^V^r • O^ttN. (41) 
The coupling constant (? p7r7r is determined from the decay width r(p — > nn), which reads 

n P ^^) = ^{M 2 p -AMlf\ (42) 

Using T(p° -> 7T+7T-) = 150.7 MeV [42], we obtain 

g p7T7T = 6.04. (43) 
Then the 27r-exchange transition amplitude with intermediate nN channel reads 
M? N = V»( q ',q)B»(k, q ') 



1 M N 1 eg pn7T g nNN , 



W -pb + p'T(pb -p-q') 



{2*y E N {p B )2E„{<f) 4M= 

" (PB -rr- M u PB -Jr- Mi mr °b ) - (44) 



where 



r = yt-j B )WB-M N )W B -t>). (45) 

The loop integration must be regularized by introducing form factors. We include the form 
factors for each vertices. In addition, we also introduce the form factor to take into account 
the off-shell-ness of the intermediate states, 

Thus the final form of the form factor is 

(47) 
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where 

A 2 — M 2 A 2 - /If 2 

TP (+\ — —1 TP (+\ — kNN 1VI n I aq\ 

and ti = (p B — p) 2 and £2 = (Pb — v') 2 ■ Here the inclusion of F pn7T (ti) implies the vector 
meson dominance assumption. The cutoff parameters will be discussed in Sect. III. 

We now comment on the loop calculation described above. We do not consider the 
crossed diagrams of Fig. 2, since such diagrams include three-particle intermediate states 
and hence are of higher-order effects which are neglected in this exploratory study. However 
by neglecting the crossed diagrams, the resulting amplitude does not satisfy gauge invariance. 
In this study, therefore, we restore gauge invariance of the amplitude (44) by projecting out 
the gauge non-invariant terms as [37] 

MT -> V^'M^V"'", (49) 
where the projection operator reads 

V V = gP- V _ (50) 



III. CROSS SECTIONS AND POLARIZATION ASYMMETRIES 

In this work we first re-examine the commonly employed a exchange by considering 
model (A) which includes the Pomeron, a, ir, i] exchanges, and the s- and -u-channel nucleon 
terms. We then explore model (B) which is constructed by replacing the a exchange in 
model (A) by the f 2 and 2tt exchanges. We also add the o exchange to model (B) as 
a correlated 2n exchange with the couplings determined by p — * 7r°7r°7 decay and Bonn 
potential. All parameters of the models are explained in Section II. In particular, the a- 
exchange parameters in model (A) are given in Eq. (9), which are chosen to reproduce the 
p photoproduction data. 

For model (B), we use the f 2 couplings as (see Appendix) 

G} nn /Atc = 2.2, F fNN = 0, G fVV = 5.76, (51) 

with the relation (29). The recently estimated T(p — > 0-7) [26] is used to constrain g pal as 

g P ^ = 0.25. (52) 

The other parameters for the a exchange are the same as given in Eq. (9). The only 
unspecified parameters are the cutoff parameters A/a^ and A/^7 for the f 2 exchange and 
the cutoff parameters of Eq. (47) for regularizing the loop integrations. The parameter A e 
for all loop integrations is fixed to be 0.5 GeV which is identical to the value used in our 
previous investigation [18] of the one-loop corrections on uo photoproduction. The other 
cutoffs including A pwn and A nNN in the loop calculation are chosen to be 0.6 GeV. The 
other two parameters of model (B) are adjusted to fit the cross section data and are found 
to be 

A fNN = A fVl = 1.4 GeV. (53) 
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FIG. 3: Differential cross sections of model (A) at £ 7 = (a) 2.8, (b) 3.28, (c) 3.55, and (d) 3.82 
GeV. The dot-dashed lines are from a exchange and the dashed lines are without a exchange. The 
solid lines are the full calculation. Experimental data are from Ref. [52] (open squares) and Ref. [3] 
(filled circles). 



This is a unsatisfactory aspect of this work, but it is unavoidable in any phenomenological 
approach. Future theoretical calculations of form factors are therefore highly desirable. 

The differential cross sections for 7p — > pp calculated from model (A) are compared with 
the SLAC data [52] and the recent CLAS data [3] in Fig. 3. We see that the full calculations 
(solid curves) are dominated by the a exchange contributions (dot-dashed curves). The 
contributions from the other exchange mechanisms (dashed curves) become comparable 
only in the very forward and backward angles. This is mainly due to the fact that the 
Pomeron exchange [Fig. 1(a)] is forward peaked and the w-channel nucleon term [Fig. 1(d)] 
is backward peaked. It is clear that the data can only be qualitatively reproduced by model 
(A). The main difficulty is in reproducing the data in the large \t\ (larger than about 3 
GeV 2 ) region. No improvement can be found by varying the cutoff parameters of various 
form factors of model (A). This implies the role of other production mechanisms in this 
region. 

The differential cross sections calculated from model (B) are shown in Fig. 4. The solid 
curves are the best fits to the data we could obtain by choosing the cutoff parameters given 
in Eq. (51) for the f'2 exchange. In the same figures, we also show the contributions from 
the f'2 exchange (dot-dashed curves), 2ir and a exchanges (dashed curves), and the rest of 
the production mechanisms (dotted curves). It is interesting to note that the ji exchange 
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FIG. 4: Differential cross sections of model (B) at E 1 = (a) 2.8, (b) 3.28, (c) 3.55, and (d) 3.82 
GeV. The dot-dashed lines are from ji exchange, the dashed lines are from 2ir and a exchanges, 
and the dotted lines are from the other processes, i.e., without /2, 2ir, and a exchanges. The solid 
lines are the full calculation. Experimental data are from Ref. [52] (open squares) and Ref. [3] 
(filled circles). 

in model (B) (dot-dashed curves in Fig. 4) drops faster than the a exchange in model (A) 
(dot-ashed curves in Fig. 3) as t increases. On the other hand, the 2n and a exchanges 
(dashed curves in Fig. 4) give a nontrivial contribution in large \t\ region at lower energies 
but are suppressed as the energy increases. Therefore such effects are expected to be seen 
at energies very close to the threshold. As expected, the a meson exchange contribution is 
much suppressed than in model (A). 

Thus we find that model (B) is comparable to model (A) that is the commonly used 
a exchange model in fitting the differential cross section data of SLAC and TJNAF. In 
particular, the data at small \t\ (< 2 GeV 2 ) can be equally well described by both models, 
as more clearly shown in Fig. 5, where the full calculations of two models are compared. 
On the other hand, both models cannot fit the data at large \t\ (> 2 GeV 2 ). But this is 
expected since we have not included N* and A* excitation mechanisms which were found 
[14] to give significant contributions to uj photoproduction at large \t\. However, we will not 
address this rather non-trivial issue here. The main difficulty here is that most of the reso- 
nance parameters associated with isospin T = 3/2 A* resonances, which do not contribute 
to uj photoproduction, are not determined by Particle Data Group or well-constrained by 
theoretical models. Before we use our model to determine a large number of resonance 
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FIG. 5: Differential cross sections of model (A) and (B) at E y = (1) 2.8, (b) 3.28, (c) 3.55, and 
(d) 3.82 GeV. The dashed lines are the results of model (A) and the solid lines are those of model 
(B). Experimental data are from Refs. [3, 52]. 

parameters by fitting the existing limited data, it would be more desirable to further test 
and improve the nonresonant amplitudes such as including more complete calculations of 
27r-exchanges. Hence, in this paper, we focus on exploring which experimental observables 
are useful for distinguishing more clearly the Model (B) from Model (A) in the small \t\ 
(< 2 GeV 2 ) region where both models can describe the differential cross section data to a 
large extent and the N* and A* effects are expected to be not important. Experimental 
verifications of our prediction in this limited t region will be useful for understanding the 
non-resonant amplitudes of p photoproduction at low energies. 

We have explored the consequences of the constructed models (A) and (B) in predicting 
the spin asymmetries, which are defined, e.g., in Ref. [37]. The results for the single spin 
asymmetries are shown in Fig. 6 for E y = 3.55 GeV. Clearly the single spin asymmetries 
including the target asymmetry T y , the recoiled proton asymmetry P y , and the tensor asym- 
metry V xxyy of the produced p meson would be useful to distinguish the two models and 
could be measured at the current experimental facilities. Of course our predictions are valid 
mainly in the small t region since the A^* and A* excitations [14] or G-pole contributions 
[22], which are expected to be important at large t, are not included in this calculation. 

Our predictions on the beam-target and beam-recoil double asymmetries [37] are given in 
Fig. 7. Here again we can find significant differences between the two models in the region 
of small \t\. Experimental tests of our predictions given in Figs. 6 and 7, therefore, will be 
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FIG. 6: Single spin asymmetries of model (A) and (B) at E 1 = 3.55 GeV. Notations are the same 
as in Fig. 5. The definitions of the spin asymmetries are from Ref. [37]. 
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FIG. 7: Double spin asymmetries CfJ, Cf?, Cf£, and Cf^ of model (A) and (B) at £ 7 = 3.55 
GeV. Notations are the same as in Fig. 5. 



useful in understanding the non-resonant mechanisms of p photoproduction. 

Since both the a and f 2 exchanges are natural parity exchanges, it would be difficult to 
test them using parity asymmetry or photon asymmetry that can be measured from the 
decay distribution of the p meson produced by polarized photon beam. For completeness, 
we give the predictions of the two models on these asymmetries in Fig. 8. As expected, 
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FIG. 8: Spin asymmetries P a and of model (A) and (B) at E 1 = 3.55 GeV. Notations are the 
same as in Fig. 5. 

it is very hard to distinguish the two models in the forward scattering angles with these 
asymmetries. 

IV. SUMMARY AND DISCUSSION 

In this paper we have re-examined the a-exchange and / 2 -exchange mechanisms of p 
photoproduction reactions. It is found that the commonly employed cr-exchange amplitude 
is weakened greatly if the a coupling constants are evaluated by using the recent information 
about the p — > 7r°7r°7 decay and the aNN coupling constant of Bonn potential. This has 
led us to introduce the un-correlated 2n exchange amplitude with nN intermediate state. 
This leading-order 2-n exchange amplitude can be calculated realistically using the coupling 
constants determined from the study of pion photoproduction and the empirical width of 

P — > 7T7T. 

In the investigation of / 2 -exchange mechanism, we evaluate its amplitude using an ef- 
fective Lagrangian which is constructed from the tensor structure of the / 2 meson. Phe- 
nomenological information together with tensor meson dominance and vector meson domi- 
nance assumptions are used to estimate the / 2 coupling constants. This approach, which is 
more consistent with the conventional meson-exchange models, is rather different from the 
/ 2 -exchange model of Laget [10], where the f 2 interaction structure was borrowed from that 
of Pomeron exchange assuming Pomeron-/ proportionality, i.e., / 2 -photon analogy. 

In comparing with the existing differential cross section data, we find that a model with 
the constructed 2n, a, and / 2 exchanges is comparable to the commonly used a-exchange 
model in which the a coupling parameters are simply adjusted to fit data. Both models 
can describe the data equally well in the small \t\ (< 2 GeV 2 ) region, but fail at large 
\t\. We suggest that experimental verifications of the predicted single and double spin 
asymmetries in the small \t \ region will be useful for distinguishing two models and improving 
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our understanding of the non-resonant amplitude of p photoproduction. 

Finally, we would like to emphasize that the present investigation is just a very first 
step toward obtaining a complete dynamical exchange model of p photoproduction at low 
energies. The next steps are to examine the additional 27r-exchange mechanisms due to, for 
example, uN and 7rA intermediate states and the crossed diagrams of Fig. 2. The effects 
due to N* and A* effects must be included for a realistic understanding of the interplay 
between the non-resonant and resonant amplitudes. Theoretical predictions of the resonance 
parameters associated with A* resonance states will be highly desirable for making progress 
in this direction. 
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APPENDIX: TENSOR MESON DOMINANCE AND / 2 -HADRON INTERAC- 
TIONS 

The free Lagrangian and the propagator of the tensor meson were studied in Refs. [53-57] . 
The propagator of the tensor meson which has momentum p reads 

Q„ = 1 ppu;pa (A.l) 

p 2 -M 2 f +ie ' y J 



where Mf is the tensor meson mass and 



with 



p„ = 1 (gwgua + _ l -g^gP^ ( A .2) 



9pu = ~9pu + Vj ^- (A.3) 



1. /27T7T Coupling 

The effective Lagrangian for f 2 irir interaction reads [45] 

Cf ^ = - Z -^d^-d^r, (A.4) 
M f 

where is the / 2 meson field. This gives the vertex function as 

Vfr = -^(Pa+Pc)M+PcU^(X f ), (A.5) 
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where p a and p c are the incoming and outgoing pion momentum, respectively. The minus 
sign in the Lagrangian (A. 4) is to be consistent with the tensor meson dominance [58]. The 
Lagrangian (A. 4) gives the f 2 — > nn decay width as 

q2 ( m 2 V /2 

r ^^ = -^ M '[ 1 - 4 w,) ■ (A6) 

Using the experimental data, T(f 2 — > 7r7r) cxpt 156.9 MeV [42], we obtain 

G 2 

w 2.64, (A.7) 

47T 

which gives Gf^ w w 5.76. 



2. Tensor meson dominance 

The tensor meson dominance (TMD) is an assumption of meson pole dominance for 
matrix elements of the energy momentum tensor just as the vector meson dominance (VMD) 
is a pole dominance of the electromagnetic current. By using TMD, one can determine the 
universal coupling constant of the f 2 meson from its decay into two pions, which can then 
be used to determine the f 2 NN and f 2 VV couplings. When combined with VMD, this also 
allows us to estimate the ^77 and /V7 vertices. It is interesting to note that the TMD 
underestimates the empirical f 2 NN coupling while it overestimates the f 2 — > 77 decay 
width. But it shows that the f 2 couplings with hadrons and photon can be understood by 
TMD and VMD at least qualitatively. Here, for completeness, we briefly review the method 
of Refs. [46, 51] to illustrate how to use TMD to get the / 2 -hadron couplings. 

Let us first apply TMD to spinless particles [46, 59]. The energy-momentum tensor 
between spinless particles can be written as 

(pl^(O)b') = Fi(A 2 )E jU E i , + F 2 (A 2 )(A M A„ - g,„A 2 ), (A.8) 

with S M = (p + p')^ and A M — (p — p')^. Then with the covariant normalization one has 

(p\ J e 00 (x)d 3 x\p) = EN P , (A.9) 

where N p is the normalization constant. By comparing with Eq. (A.8), one can find 

*i(0) = \ (A.10) 



Now we define the effective couplings for tensor mesons as 



</IM°)l°> = 9 f Mfe, u , (p\f\p') = - C ^E"2^, (A.ll) 
where the latter equation is consistent with Eq. (A. 5). The pole dominance gives 
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ay/3 Gfpp 



Mf A 2 - Mj 



f 

Til 



which leads to 



Thus we have 



_ g f MjG fpp 1 1AA 2 

" A 2 - MJ 3^ L + 3 M 2 L J ' (A - 12) 

/ 1 

It should be noted that the sum of Eq. (A. 14) contains tensor meson nonet, i.e., /2(1270) 
and /2(1525). But in the case of the fi^is coupling, if we assume the ideal mixing between 
the / 2 (1270) and the / 2 (1525), the /a(1525) decouples by the OZI rule. Therefore we obtain 
Gfi nw ~ 0, and the universal coupling constant is determined as 

^ = -J— » 0.087, (A.15) 
ZLr f™ 

using the value of Eq. (A. 7). 

With the universal coupling constant determined above, one can now use it to estimate 
the f 2 NN coupling. For this purpose, we apply TMD to spin- 1/2 baryon state. The energy- 
momentum tensor of the spin-1/2 baryons can be written as 

(p|M<W> = u(p){^(7^ + 7*E,.)*i(A 2 ) 

+ lM7 F2(A2) + (A " A " " ^ A2 ) F 3(A 2 )} u{p'). (A.16) 
With the covariant normalization, the conditions 
(p| J 0oo(x)d 3 x\p) = EN p , 

(p,p = 0,s 3 = +\\${x 1 9o 2 {x) - x 2 i(x)}d 3 x\p, p = 0, s 3 = +|) = ±N P , 

(A.17) 

give 

F!(0) = 1, F 2 (0) = 0. (A.18) 

Now using the form for f 2 NN coupling in Eq. (26), assuming the pole dominance gives 
the following relations: 

M f M r 
-1 = Ag J G fNN —+lg J ,G f . NN — , 

M f M f , 
= 4g f Ff NN ^ + 4g f ,Ff, NNJ J-. (A.19) 
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Again by assuming the decoupling of the f 2 from the nucleon coupling, we can have [46] 



_ 1 M p _ G f7T7T M p 

G -f NN ~ 4f f W f ~ ^TW f ~ 2A2 ' 

F fNN = 0. (A.20) 

This gives G 2 NN /4ir ps 0.38 as shown in Table I, which is smaller than the values estimated 
by nN dispersion relations by an order of magnitude. It should also be noted that the 
values estimated by nN dispersion relations may be affected by the inclusion of other meson 
exchanges. More rigorous study in this direction is, therefore, highly desirable. 



3. fiVV coupling 

Before we discuss /277 and ^2^7 couplings, we first apply TMD to f^V coupling, where 
V stands for vector mesons. The energy-momentum tensor between identical vector mesons 
contains six independent matrix elements [51], 

(V\e^\V'} = g 1 {A 2 ){e-e')T 1 ^ v + g 2 {A 2 ){e-T 1 ){e' • E)E M E, 

+ Q 3 (A 2 ) {(e • T)e'^ v + (e • + (e' • £)e M E„ + (e' • E)e v E ll } 

+ £ 4 (A 2 ) {(e • A) e ;A, + (e • A)e^A M + (e' • A)e (ti A l/ + (e' • A)e,A M 

-2(e • A)(e' • A)^ - A 2 (e fl e' l/ + &„)} 
+ e 5 (A 2 )( e .e')(AA-AV) 

+ £ 6 (A 2 )(e • S)(e' • E)(A /Lt A l/ - ^,A 2 ), (A.21) 

where E M = (p + A M = (p — p')^ and e, e' are the polarization vectors of V and V, 
respectively. Then the conditions like Eq. (A. 17) give 

&(0) = ~, &(0) = \ (A.22) 

In the pole model, the form factors £i(A 2 ) . . . £ 4 (A 2 ) are dominated by tensor meson 
poles. Because of the symmetry property of the tensor meson, we have generally four /2VV 
coupling vertices: 

(V\f\V) = ^(e-6')(E /1 E^) + ||(6-S)(6'.E)(E,E^) 
C 

^3 f I _ v-iN _/ v-i , / _ VI \ _/ V , I J V\ _ v-i , / J Vl\ . X-l 1 fflV 



+ ^(-A 2 )(6,6'+6^)/^, (A.23) 

while we have used A^f^ = f£ = in writing the G 4 term. For our later use, an effective 
vertex H(A 2 ,p 2 ,p' 2 ) is introduced to replace — jjjA 2 as [51] 

H(A 2 ,p 2 , P ' 2 ) = 2L {-A 2 + a( P 2 + P ' 2 - 2M 2 )} . (A.24) 
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FIG. 9: f 2 — > 77 decay in vector meson dominance. 



Now we use the pole dominance again using Eq. (A. 11) to find 



3fMjG 



e ' (A2) = iSf < A - 25 > 

which leads to 

l - = g f G u G 3 = -d, (A.26) 

combined with Eq. (A. 22). Therefore, with Eq. (A. 15) we get 

G l = -G 3 = G f7r7T « 5.76. (A.27) 

The above relation should hold for f 2 pp and f 2 uu. The SU(3) symmetry and the ideal 
mixing give 

G 1 (&<M>) = V2G 1 (f 2 pp), G 1 (&<f>u) = G 1 (f 2 <M>)=0. (A.28) 

Note that two couplings G\ and G3 are determined by TMD but G 2 and G4 cannot be 
estimated without further assumptions. 

4- /277 and f 2 V~y couplings 

The remaining two couplings G 2 and G 4 of Eq. (A. 23) are estimated by using VMD and 
gauge invariance. We consider f 2 — > 77 using VMD as illustrated in Fig. 9. 
By using e • k — e' • k' — and VMD, we have 



2 

e 



<7(*)7(W> = ( p_ M ^-M*) {t^^ 

- ||(e • k'w ■ k)(k - k'uk - k\r 

+ ^[-(e-k%(k-k%-(e.k'K(k-ki)„ 

+ (e' ■ k)e^{k - k') v + (e' • k)e v (k - fc%] r 
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with 



where we have introduced the notation 

/ M 2 \ 2 / /If 2 \ 2 

Gi = \-JL) G f / P + G f r, (A.30) 

M 2 

(QbTW = -r^{v). (A.3i) 

jv 

Because of isospin, there is no mixing between the intermediate p and uo mesons. By 
looking at the amplitude (A. 29), however, one can find that it is not gauge invariant, i.e., it 
does not vanish when replacing e M by k^. This gives a constraint on the couplings. The most 
general form for f 2 j~f satisfying gauge invariance has two independent couplings as [51] 

-,2 f 



(l(kHk')\f) = -^{A [(e ■ e')(k ■ k') - (e ■ k')(e> ■ k)} (k - k\{k - k%f 

+ B[(e ■ e')(k - fc%(fc - k') v + e'^k - k') v {e ■ k') + e' u (k - fc%(e • k!) 
- e,{k - k') v {e' ■ k) - e v {k - k')„(e' ■ k) 

-2(k-k')(e^ + e'^)]r}, (A.32) 



which then gives 



G2_ 



(k-k')A + B, 

A, 

-B, 



[-Mj + a(k 2 + k' 2 - 2M 2 )] = -2(k ■ k')B = (k 2 + k' 2 - M 2 )B. (A.33) 

Solving this system at k 2 = k' 2 = and G± = —G3 gives 

A = G 2 = 0. (A.34) 

Since gauge invariance applies to jsoscalar and isovector photons separately, we get G2 = 
for V = p, uj. Still we do not fix G 4 and a, but have a constraint, 

G 4 (M 2 f + 2aM v ) = G\M 2 . (A.35) 

To complete the model, let us finally consider /V7 vertex. Here again, we use the VMD 
as in Fig. 10. The gauge invariance of the vertex at k 2 = M v and k' 2 = leads to 

G 4 {M 2 f + aM 2 ) = G 1 (M 2 f - My). (A.36) 

Then solving the coupled equations (A.35) and (A.36) gives [51] 

. Ml 2 Mr Mi 
G A = G l f -— 2 v , a= —„ f. . (A.37) 
M 2 Mf - 2.1/ 2 v ' 
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FIG. 10: /2 — > decay in vector meson dominance. 

Thus we have determined all couplings of Eq. (A. 23) with the relation (A. 30). 

The above procedure shows that the /277 and /2^7 vertices can be written with two 
form factors because of gauge invariance, which read 

< 7 (*)7(W2> = ^' x F u A{i; u (k,k'), 

<7(*)V(A0|/ 2 ) = ^ K e ,x rA{ v ^(k,k'), (A.38) 

where 

AlHAk, k') = ^ [g KX (k ■ k') - k'M (k - k%(k - k') v 

+ 9f 71 [9K\(k - k')^{k - k') v + gx^K( k ~ k ')v + 9xuK( k ~ k ')v 
9Kiik\{k k) u Qnvkxik k )^ 
- 2k ■ k'(g Kf ,g X u + SWfiv)]- (A.39) 

The vertex function A^^ v (k, k') can be obtained from A^ 1 by replacing j ^ n and g^ by 
g f vy and 0/y 7 » respectively. 

With Eqs. (A.38) and (A.39), we can obtain the f 2 — ► 77 decay width as 3 

r (/ 2 -77) = ^(^// 2 77 + ^ 77 ). (A.40) 
Then TMD and VMD give [51] 

2/I 1 



//T7 = °' ^/77 = e \Jj + Jl) GfVV - (A ' 41) 

The vector meson decay constants are f p = 5.33, / w = 15.2, and = 13.4. By noting that 
TMD gives G/yy = Gf n7n we get 

r(/ 2 -> 77) « 8.8 keV, (A.42) 



Here we do not agree with the decay width formula of Ref . [60] . 



24 



while its experimental value is T(/ 2 — > 77) exp t. = 2.6 ± 0.24 keV. Thus we can find that this 
procedure overestimates the experimental value by a factor of 3 ~ 4. 

The decay width of fi — > can be computed using Eqs. (A. 38) and (A. 39) as 

M fn ^3 f 1 ,2,, ^4 * , « „ ,x(l-x) 2 



r(/ 2 -^7) = ^(i-^) 3 {^IV 7 l 2 (i-^) 4 -(W^ 7 + //V 7 W 



12 



2 



2 6 

where x = M^/M 2 f . TMD combined with VMD gives [51] 



//v 7 = 0, 5 /y7 = (A.44) 



This leads to 

9fPi = ft 

and [58] 



r(/ 2 - P7)/r(/ 2 - c 7 ) = = =g = 8.14 ± 1.2 (A.45) 



r(/ 2 - P7)/r(/ 2 - 77) = 2geZ(l - x) 3 (l + f + yj = 155. (A.46) 

Those quantities are not measured yet. Therefore, measuring those quantities will be very 
useful to test TMD and VMD. 
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